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Abstract 

We aim to construct from first principles a perturbative framework for studying nonequi- 
librium quantum-field systems that include massless Dirac fermions. The system of our 
concern is quasiuniform system near equilibrium or nonequilibrium quasistationary system. 
We employ the closed-time-path formalism and use the so-called gradient approximation. 
Essentially no further approximation is introduced. We construct a fermion propagator, 
with which a well-defined perturbative framework is formulated. In the course of construc- 
tion of the framework, we obtain the generalized Boltzmann equation (GBE) that describes 
the evolution of the number- density functions of (anti)fermionic quasiparticles. 

11.10. Wx, 12.38.Mh, 12.38.Bx 



I. INTRODUCTION 

Ultrarelativistic heavy-ion-collision experiments at the BNL Relativistic Heavy Ion Col- 
lider (RHIC) have begun and will soon start at the CERN Large Hadron Collider (LHC) 

*Electronic address: niegawa@sci.osaka-cu.ac.jp 
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in an anticipation of producing a quark-gluon plasma (QGP) [see, e.g., The QGP to 

be produced is an expanding nonequilibrium system. Studies of the QGP as such have just 
begun. 

In previous papers, a perturbative framework has been formulated from first principles 
for dealing with out-of-equilibrium complex-scalar field system [^], O(N) linear-sigma system 
@], and the system that includes gauge bosons ||. In this paper, we take up the out-of- 
equilibrium quantum-field theories that include Dirac fermions. Keeping in mind that the 
masses of light quarks may be ignored in QGP, we treat massless fermions. [Generalization 
to massive ones is straightforward.] Through similar procedure as in we construct 

the fermion propagator and, thereby, frame a perturbation theory. Essentially, only ap- 
proximation we employ is the so-called gradient approximation (see below). We use the 
closed-time-path (CTP) formalism 0-§J of nonequilibrium statistical quantum-field theory. 

Throughout this paper, we are interested in quasiuniform systems near equilibrium or 
nonequilibrium quasistationary systems. Such systems are characterized (cf., e.g., ||) by 
two different spacetime scales: microscopic or quantum-field-theoretical and macroscopic 
or statistical. The first scale, the microscopic-correlation scale, characterizes the range of 
radiative correction to reactions taking place in the system while the second scale measures 
the relaxation of the system. For a weak-coupling theory, in which we are interested in 
this paper, the former scale is much smaller than the latter scale. A well-known intuitive 
picture || for dealing with such systems is to separate spacetime into many "cells" whose 
characteristic size, L M (/i = 0, ...,3), is in between the microscopic and macroscopic scales. 
It is assumed that the correlation between different cells is small, so that microscopic or 
elementary reactions can be regarded, to a good approximation, as taking place in a single 
cell. On the other hand, in a single cell, prominent relaxation phenomena do not take place. 
This intuitive picture may be implemented as follows. Let S(x, y) be a propagator. For 
a system of our concern, S(x,y), with x M — y M [|x M — y^\ ^5 L M ] fixed, does not change 
appreciably within a single cell. Thus, the mid-point X M = (x M + y M )/2 may be used as a 
label for the spacetime cells and is called the macroscopic spacetime coordinates or a slow 
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variable fLOfl . On the other hand, relative spacetime coordinates x M — y^, which is called 
a fast variable, are responsible for describing microscopic reactions taking place in a single 
spacetime cell. We introduce a Wigner transformation (Fourier transformation with respect 
to the relative coordinates x — y with (x + y)/2 held fixed): 

S{x,y) = J ^e- iP <^S(X;P) (1.1) 

(with P M = (p°,p)). The Wigner function S(X;P) depends on X only weakly. The self- 
energy part E(x, y) enjoys a similar property. We shall freely use S(x, y) or S(X; P) [S(x, y) 
or E(X; P)], which we simply write S [E] whenever obvious from the context. 
In the following, there often appears a "product" of functions: 

[F-G}(x,y) = Jd 4 zF(x,z)G(z,y), (1.2) 

For the Wigner transform of the function [F ■ G](x,y), we use the gradient approximation, 

[F ■ G](X; P) = Jd\x- y) e lP ^[F ■ G](x, y) 
~ F(X; P)G(X; P) 

- 1 -{F(X;P), G(X;P)} PB , (1.3) 

where X = (x + y)/2 and 

dF(X; P) dG(X; P) 



{F, G} PB 



dX» dP^ 
dF(X; P) dG(X; P) 
dP^ dX^ 



;i.4) 



We refer to the first term on the right-hand side (RHS) of Eq. ( |1.3| ) as the leading part 
(term) and to the second term, which includes derivative dx, as the gradient or nonleading 
part (term). Crudely speaking, the leading term is the quantity in itself in the spacetime 
cell labeled by X, while the gradient term represents the effect arising from the (weak) 
correlation between the spacetime cell labeled by (x + z)/2 (see Eq. (|1.2p ) and the one 
labeled by (z + y)/2. 
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The perturbative framework to be constructed accompanies the generalized Boltzmann 
equation (GBE) for the number density of quasiparticles. The framework allows us to 



compute any reaction rate by using the reaction-rate formula deduced in [[LI]]. Substituting 
the computed net production rates of quasiparticles into the GBE, one can determine the 
number densities as functions of spacetime coordinates X, which describes the evolution of 
the system. 

For the sake of concrete presentation, we take up QCD. However, the procedure has 
little to do with QCD and then the framework to be constructed may be used for any theory 
that includes massless Dirac fermion(s) with almost no modification. The plan of the paper 
is as follows: In Sec. II, making use of the free quark (fermion) fields in vacuum theory 
and the quark-distribution function at the initial time (X° = — oo), we construct the bare 
quark propagator, with which the "bare-iV scheme" may be constructed. The perturbative 
calculation on the basis of this scheme yields 0] divergence due to pinch singularities. Then, 
in Sec. Ill, we set up the basis for formulating the "physical-iV scheme" by introducing a 
new quark-distribution function and a new "free" quark fields. In Sec. IV, we make up 
the self-energy-part resummed quark propagator. In Sec. V, imposing the condition on the 
quark-distribution function that there do not appear large contributions, which stems from 
the above-mentioned pinch singularities, we construct a "healthy" perturbative framework. 
It is shown that, on the energy-shell, the condition turns out to be the generalized Boltzmann 
equation. In Sec. VI, we frame a concrete perturbative framework. Section VII is devoted 
to summary and discussion, in which comparison with related works is made. Concrete 
derivation of various formula used in the text is made in Appendices. 



II. CLOSED-TIME-PATH FORMALISM 
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A. Preliminary 



We follow the procedure as in ||. The CTP formalism is reduced to a two-component 
formalism 0. Every field, say 0, is doubled, — > (0i, 2 ), and the classical action turns out 

+oo p 

dx dx[£(</) 1 (x) 1 ...)-C(<l> 2 (x) 1 ...)] (2.1) 

-oo J 

with C the Lagrangian density of the theory under consideration. For definiteness, we take 
up massless QCD. The gluon sector has already been studied in ||. In this paper, we deal 
with the quark sector, 

C = i$ a (x)0ip a (x) + (2.2) 

where "a" is a color index. Perturbation theory is formulated in terms of propagators, 
vertices, and initial correlations^. The aim of the present paper is to construct the quark 
propagator, which takes the 2x2 matrix form 0, 

S%{x,y)= . (2.3) 

^ -i(mx)W\{y))e) -i(»)).(»> ) 

Here (...) = Tr(...p) with p the density matrix at the initial time, the 'caret' denotes the 
2x2 matrix as indicated on the RHS, the suffix a {(3) denotes the component of tpj (ijjj) 
[j = 1, 2], and T (T) is the time-ordering (antitime-ordering) symbol. Noting that the state 
realized at a heavy- ion collision is color singlet, we restrict ourselves to the case where p is 
color singlet, so that S^g = 5 ab S a p. We also assume that p commutes with baryonic charge. 
[Generalization to the color- nonsinglet case is straightforward.] In the sequel, we drop the 
color index. At the end of calculation we set = ip2 and ipi = 4>2 J/J. Instructive examples 
of p, in the case of scalar field, is given in |J. 



T Given a density matrix, which characterize the system, the initial correlations are determined 
PJlO[|. In perturbation theory, the initial correlations are treated as vertexes. 
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B. Quark propagator 



We start with computing Eq. (|2.3D with ipi = ip 2 (= an d tpi = ^2 (= V>)- For ip 
and we use the plane-wave decomposition with helicity basis in vacuum theory, which, in 
standard notation, reads 



J ■ 1 (zn^zp f=± 



<i 3 p 
2^)3 

+ dj(p)v £ (p)e ip -] (p°=p), (2.4) 

where £ (= ±) denotes helicity, and b^(p) [bg(p)] is the annihilation [creation] operator of a 
quark [an antiquark] with helicity £. Substituting Eq. ( [2.4|) and its adjoint into Eq. (p7^), 
we compute the Wigner transform, S(X;P), of S(x,y) [cf. Eq. (|I~T|)1 to the gradient ap- 
proximation. There emerge four types of terms, which contain (0^0^), (b^d^), and 
(b^d^), respectively. For an equilibrium system, only and (d±d±) terms survive, and 
(b+b + ) = (b_b-) and {d) + d + ) = (d_d-) hold. With this observation in mind, we see that, for 
the systems of our interest, the leading term of S(X; P), which we write S Q (X; P), contains 
the "even" combinations X^=±(^£) an d J2^=±(d\d^) . It turns out that the gradient term 
Si(X;P) contains also X^=±(^^) and ^=±(^^) together with spacetime (X) derivative. 
We assume that the term 5*2 (X; P) that contains the "odd" combinations X^=±£(^£) an d 
J2^=±^,(d\d^) and the term S3(X;P) that contains (b±b T ) and (d±d T ) are nonleading and 
are, at most, of the same order of magnitude as the gradient term S\(X; P). We also assume 
that the terms with (b^d^) and with (b^d^) may be ignored. It is a straightforward task to 
obtain the form of S(X; P) for the case where some or all of the above assumptions do not 
hold. 

Thus, S(X; P) consists of four terms, the leading term So, the gradient term Si, and the 
nonleading terms S2 + S3, 

S(X; P) = S (X; P) + SxpT; P) + S 2 (X; P) + S 3 (X; P). (2.5) 

The form for S 2 and S 3 are displayed in Appendix A. The terms S and S\ reads 



Si(X;P) 



A R (P) 

^A R (P)-A A (P) -A A (P) 
Ofiup _|_ e 3^up df(X; P) 



Here e 0123 = 1, 75 = i^y ^ 1 ^ 2 ^ 3 , and 



f(X;P)J>[A R (P)-A A (P)]A, 
[A R (P)-A A (P)]A + . 



(2.6) 
(2.7) 



A, 



±1 1 
1 



(2.8) 



f(X;P) = 6(-p ) + e(p )N(X;P), 

N(X;P) = i [%°) {iV ++ (X; \p°\, p) + i\L_(X; |p°|,p)} 

+0(-p°) {iV ++ (X; -p) + i\L_(X; |p°|, -p)} 
= %>(X; \p°\, p) + 0(-p°)n(X; -p), 



(2.9) 



(2.10) 



with 



%(X;bo|,p) = /rf 3 ge-^ x Tr[^( p - q /2)6 c (p + q/2)p] (£,C = +,-)• (2.H) 

Here |po| = P, P = p/p, and Q = (?°,q) with g = q • p/p- ^'s are antiquark counterpart 
of iV's. For the system of our concern, Tr[...] is different from from zero |f7[,|l0f only for g 
small compared with p. It can readily be seen from Eq. ( |2.11|) that when p is translationally 
invariant N^(X; |p°|,p) is independent of X. In passing, for a charge-conjugation-invariant 
system, N(X; P) = N(X; -P) holds. 
From Eq. (gXtfl follows 

%°)P-5xiV 5C (X; |/|,p) 

= 9(-p°)P ■ d x N (( (X; \p% -p) = (bo| = p). 

(2.12) 



Using this in Eqs. ( |2.9|) and Q2.10| ), we have 
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P-d x f(X;P) = (\p \=p). (2.13) 

As is obvious from the construction, or as can be directly shown, to the gradient approxi- 
mation, 

iT 3 x S(x,y) = -iS(x,y) $ y f 3 = 16 A (x - y), 

where 1 is the 2x2 unit matrix and % is the third Pauli matrix. Two equations in Eq. Q2.12 ) 
are "free Boltzmann equations." One can construct a perturbation theory in a similar 
manner as in ||. We call the perturbation theory thus constructed the bare-iV scheme, 
since N obeys the "free Boltzmann equation." Perturbative computation within this scheme 
yields divergences due to pinch singularities. In J|, how to deal with these divergences is 
discussed and shown is that the bare-iV scheme is equivalent to the physical- N scheme, to 
which we now turn. 



III. CONSTRUCTION OF THE PHYSICAL- N SCHEME 

Following the procedure as in M , we construct a scheme in terms of the number density 
that is as close as possible to the physical number density. To this end, first of all, we 
abandon the "free Boltzmann equation" (|2.13|) . This means that / here and in the sequel 
is different from / (= f B ) in the last section. Then S(X; P; f) ^ S(X] P; which means 
that ip, and ip in the present (physical- A) scheme differs from the free field in vacuum theory, 
Eq. ( |2.4|) , as employed in the bare-A scheme. [See also the comment below in conjunction 
with Eq. (^5)-] Specification of / and then also of ip and are postponed until Sec. VI. 

Now, S(x,y) is not an inverse of if 3 [cf. Eq. (JO) with Eq. (p^)]- Straightforward 
computation using the inverse Wigner transform of S(X; P), Eq. ( |2.5| ), yields, to the gradient 
approximation, 

if 3 0J(x, y ) = i5\x-y)- % -r 3 A + J ^ e -^-(-») 
x^^H^R-^A)P-d x mP). 

(3.1) 
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Our procedure of constructing a consistent scheme is as follows: We further modify if) and 
■0 by adding a suitable S a dd to S on the left-hand side (LHS) of Eq. (|3.1[). The conditions 
for Sadd to be satisfied are 

• Sadd vanishes in the bare-iV scheme in Sec. II. 



To the gradient approximation, Eq. (|3.1|) turns out to 



(tfJl-L c )-(s + S add )=ll. 



(3.2) 



Here we have used the short-hand notation ( |1.2| ). In Eq. Q3.2| ), '1' is the matrix in 
spacetime whose (x, y)-component is <5 4 (x — y), and L c is some (4 x 4) <g> (2 x 2) matrix 
function. 

It is straightforward to obtain the form of the required S a dd'- 

7 ° + 7 3 P 



S add (X; .P) = iA, 



/\ 2 4- 



p° + p 3 P 2 



xP-d x f(X;P), 



(P/P 2 the principal part of 1/(P 2 ± *0 + )) from which we obtain for L c in Eq. 

L c (x,y) = L C A. = iA_ [ i^ e - iP -^)^±4 
V ,y; J (2tt) 4 p°+p 3 

xP-d x f(X;P), 



(3.3) 



with A_ as in Eq. (|2.8|). In obtaining Eq. ( |3T2| ) with S^dd as in Eq. (|3 . 3|) , we have used 
L c ■ S ~ L c ■ (S + Sadd), since the difference can be ignored to the gradient approximation. 

In a similar manner, we find that (S + Sadd) ■ (—^1 $ —L c ) = 11. Thus we have found 
that the (2 x 2) <g> (4 x 4) matrix propagator (S + S a dd) is an inverse of (if 3 01 — L c ), so that 
the free action is 



d 4 x d A yi[i(x) if 3 x S 4 (x - y) - L c {x,y)AJ\ 0(y), 
i = ("01, $2) , 



(3.4) 



This is different from the free action in the bare- A scheme, / d 4 x i/iif^i/), and then the free 
fields in Eq. (ETj) is different from the free fields employed in the bare-A scheme. Since 



the term with L c (x, y) in Eq. (|3.4f) is absent in the original action, we should introduce the 
counter action to compensate it, 



A c = J d 4 xd / 'y^(x)L c (x,y)A_^(y), (3.5) 

which yields a vertex 

d 4 P _- iP , x -vM° + l 3 



rL c {x,y)A„ = -A_j—e-^-yy- 



(2tt) 4 p°+p 3 
xP-d x f(X;P) 

= -i£U( x ,y). (3.6) 

As in ||, it is straightforward task to construct a L c -resummed propagator, which we do 
not reproduce here. 

In closing this section, we emphasize that f(X; P) in the present scheme is an arbitrary 
function, provided that /(Jf°„, X; P) [= 0(-p°) + e(p°)N(X° n , X; P)], with Xf n (= -oo) the 
initial time, is a given initial data. We have introduced the counteraction A c , Eq. (|3.5|) , so 
as to remain on the original theory. Thus, it cannot be overemphasized that the schemes 
with different /'s are mutually equivalent. If we choose for /, which subjects to the "free 
Boltzmann equation" ( |2.13| ), the scheme reduces to the bare-iV scheme in the last section. 
In Sec. VI, we shall choose /, with which a well-defined perturbation theory is formulated. 
As will be shown below, it turns out that / is determined order by order in perturbation 
theory. As a natural assumption, we require that, in the limit g — > (g the QCD coupling 
constant), / — > fs- 

IV. RESUMMATION OF THE SELF-ENERGY PART 
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A. Preliminary 



As has been mentioned at the end of last section, when interactions are switched off, 
the self-energy part E vanishes and Eq. (|2.13|) holds, P ■ dxf(X] P) = 0. Thus, we suspect 
that the self-energy part causes the nontrivial evolution of /, P ■ dxf{X; P) ^ 0. Thus, E 
ties to the gradient part, S\, of S. More precisely, E is of the same order of magnitude as 
(S'o)~ 1 'S'i('S'o)~ 1 and (5'o)~ 1 S' a dd(5'o)~ 1 . Hence, in computing E in the approximation under 
consideration, it is sufficient to keep the leading part (i.e., the part with no A M -derivative). 

Since we are dealing with the massless quark, the leading part E(X; P) may be decom- 
posed as E = [a+P- - o-p + }/2 + Yll=i °al • e a (p), where 

P± = (l, ±p) (p = p/p), (4.1) 

and e*i(p) and e 2 (p) are unit vectors being orthogonal to p; p ■ e a (p) = ei(p) ■ e 2 (p) = 
(a = 1,2). For an equilibrium system, <5\ — &i = 0. Then it is natural to assume that, for 
the systems of our concern, a a (a = 1, 2) can be ignored when compared to the leading part, 



E(X; P) = \ [a + {X; P)$>_ - &~(X- P)?> + ] 



2 

Generalization to the case where this assumption does not hold is straightforward. It should 
be noted that E consists of two pieces, 

E = E loop + E (c) , 

where E loop is the contribution from loop diagrams and E^ c ) is as in Eq. (|3.6| ). It should be 
remarked that some £ loop contains internal vertex(es) iE^. [See Sec. VI below.] 

Within the gradient approximation, it is sufficient to perform a E-resummation for the 
leading part Sq. This is because the corrections to other parts due to the resummation are 
of higher order. Thus, for Si — S3 and S a dd, one can use the formulae in the bare-iV scheme 
in Sec. II [see the argument at the end of Sec. III]. In particular, for / in S\ — S3 and S a dd, 
one can use Now we introduce, as usual, the "standard form" (see Eq. (|1.2|)) p2|fl 
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B L ■ Sdiag • B R (x,y), 
Sdiag = diag(i0A R , -i$A A ), 



B, 



v 



B 



R 



J 



(l 
\ 



-f ~f 
1 1 



(4.2) 



(4.3) 



1 l-f 

where / (= f(x,y)) is the inverse Wigner transform of f(X;P) and l(x,y) 
Computing Eq. (|4.2|) to the gradient approximation, we obtain 



5 i (x - y). 



S(x,y) = S {x,y) + S^x^) 



B L ■ -S'diag • B R (x, y) + i 



d A P 

Wr 7 



-iP-(x-y) 



> x f(X;P)-^±^P-d x f(X;P) 

pU _|_ p6 



(4.4) 



It is obvious that one can freely include gradient part(s) into the "resummed part." For 
convenience, we include the gradient part S'add, Eq- (I3.3Q , and take S in Eq. (|4.4| ) as the 
"resummed part." 

It is to be noted that, from Eqs. ( pH|) and (|3.3|) follows^ 



E R i+i <% = E H i+j 

i,j=l i,j=l 



Sij + (^add, 



'.I 



0. 



(4.5) 



B. Self-energy-part resummed propagator 



A S-resummed propagator G obeys the Schwinger-Dyson (SD) equation: 



G = S + S-Z-G = S + G-Z-S. 



(4.6) 



We recall that G obeys |7j the same relation as Eq. §J), Y? itj=l {-) i+j Gij = 0. Using this 
and Eq. flj.5|) in Eq. (|4.6| ), we obtain 



^Here and in the following, (ij) element of a (2 x 2) matrix M is denoted by Mij. 
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(4.7) 



Procedure of solving Eq. ( |4.6| ) is given in Appendix B. The result for G(X; P) takes the 
form 



G{X-P) 
( 



G R (X;P) 
^ G R (X; P) - G A (X; P) -G A (X; P) 
- A + [G R (X; P) - G A (X; P)} f(X; P) 

+ A + G K (X;P). 
Gr(a) consists of the leading piece and the nonleading piece, 



G R{A) (X; P) ~ G^ A) (X; P) + G^ (A) {X- P), 



R(A)\ 



The leading piece reads 



G r}a)( X i P ) — 2 



9^A)(X;P)h + 9 R{ AM;P)^ 



9R(A) ~ o T „ =r n ± ' 
V +P + a R(A) 



(4.8) 



(4.9) 



(4.10) 
(4.11) 



where P± is as in Eq. (|4.1| ) and & R (a) = a li + cr i2(2i) • Although the nonleading piece G^j A \ 
can be ignored to the gradient approximation, we have displayed it in Appendix B. 
G K in Eq. (fU3) reads 



G K ^G$ + G$, 



G 



(o) 

A 



£ {(P T ■ d x f - ir (^) lo ° p 

+r{Rea R , f} P . B .) g R g T J T 



T=± 



(4.12) 



(4.13) 



G * = ^ £ £ [{Pr ■ dxf 



H=R,At=± 



+t{* t h, f} P .B.)(9 T n) t 
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16 p 4 



where 7^ = 7* — (7 • p)p* and 



- [1 - /] (^P + / (4) 



loop 



\ loop 



(4.14) 



(4.15) 



In obtaining Eq. ( [4.14 ), use has been made of (o'a) = a Ri which can straightforwardly be 
proved H. 



V. GENERALIZED BOLTZMANN EQUATION 



A. Energy-shell and physical number densities 



For later use, referring to Eq. ( 4.10 ) with Eq. fl4.11|) , we define the energy-shell for 
"normal modes" through 



Re 



p°TPTRea±(X;P) 



p°=±UJl(X;±p) 



0. 



(5.1) 



It is well known that, in equilibrium quark-gluon plasma, "abnormal modes" called 
plasmino appears for soft p (= O(gT)) [g the QCD coupling constant and T the temperature]. 
The energy-shell of such modes, if any, is defined through 



Re 



pO= T o;|(X;Tp) 

p = F p= F Reai{X ] P) 



P 0=T^(X;=FP) 



0. 



(5.2) 



Useful formulae that hold on the energy-shell are displayed in Appendix C. 

In order to obtain the expression for physical number densities, we start with computing 
charge and momentum densities, 
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Tr 



f{x)p 



Tr 



tp(x) r Y°il>(x) p 



Tr 



V(x)p 



— Tr 

2 



-zTr 

±Tr 
2 



7° (G 2 i(x, x) + G 12 (x, x j) p 



7°V (G 21 (x,x) + G 12 (x,x)) p 



(5.3) 



(5.4) 



We first compute the leading contribution to Tr [j°(x)p\. Substituting the leading parts 
of G 21 and of G X2 (cf. Eq. (pTg|)) and using Eq. ( ft.lUj ), we obtain 

d A P 



Tr 



f(x)p 



(g*(x;P)+g?(x;P))-c.c. 
x [2/(x;P)-l]. 



From this, we see that Tr [j°(a;)p] weakly depends on x. The narrow-width approximation, 
Ima± — > =)=i0 + yields 



Tr 



j°(x)p 



d 3 p 
(27p 



Z"(x; w?(p), p)n(z; <(p), p) - Z n _{x- u£(p), p)n(x; c^(p), p) 



+Z a _ (x; ^(p), p)n(x; ^(p), p) - Z« (x; w° (p), p)n(x; ^(p), p) + .... (5.5) 



Here n and n are as in Eq. (|2.10| ), and stands for the contribution from 2f — 1 3 —e(p°) 
[ cf. Eq. (|2.9| )1, which is the vacuum-theory contribution corrected by the medium effect. 
Z's in Eqs. (|5.5|) are the wave-function renormalization factors, Eqs. ( |C.1| ) and ( |C.4| ) in 
Appendix C. The first (last) two contributions on the RHS come from the "normal modes" 
("abnormal modes") of quasiparticles. If there are several normal and/or abnormal modes, 
summation should be taken over all modes. The factor '2' in Eqs. ( |5.5| ) comes from the spin 
degrees of freedom. 

In a similar manner, we obtain for the momentum density, 



Tr 



V{x)p 



d 3 p 
(2vr) 3 



P 



Zl{x; ul(p),p)n(x; w?(p), p) + Z%x- u£(p), p)n(x; u£(p), p) 



-Z a _{x- u£ (p), p)n{x; u a _(p), p) + Zl(x; u a Jp), p)n(x; u a Jp), p) + .... (5.6) 
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From Eqs. ( |5.5|) and (|5.6|) , we can read off that n is the number density of fermionic quasi- 
particle and n the number density of antifermionic quasiparticle. Undoing the narrow-width 
approximation yields further corrections to the physical number densities. 

Let us turn to analyze the contributions from the nonleading part of G in Eq. (JO|). 
Inspection of Eqs. Q and (^) with Eqs. Q - ([03]) and P^DQ shows that all but G%\ 
Eq. (|4.13|) yield well-defined corrections to the physical number densities due to the medium 
effect. G$ contains 

9r9a ~ 



p° =F p =F ctr 



In the narrow-width approximation Imcr± — > =F0 + , g±g± develops pinch singularity in a 
complex p°-plane. Then the contributions of G% to Eqs. ( |5.3| ) and ( |5.4| ) diverge in this 
approximation. In practice, Imcr^ (oc g 2 ) is a small quantity, so that the contribution, 
although not divergent, is large. This invalidates the perturbative scheme and a sort of 
"renormalization" is necessary for the number densities |§. This observation leads us to 
introduce the condition G^ = on the energy-shells: 



P ± ■^ x fT^(^) lOOP ±{Reat /} 



on energy shells 



= 0, (5.7) 

where {..., ...}p B is as in Eq. ( |1.4| ). This serves as determining equation for so far arbitrary 
/. (See below, for more details.) Now the above-mentioned large contributions, which 
turn out to the diverging contributions (due to the pinch singularities) in the narrow-width 
approximation, do not appear. Thus, the contributions from G$ to Eqs. ( ^3|) and to ( |5.4|) 
also yields well-defined corrections to the physical number densities. 



B. Generalized Boltzmann equation 

We are now in a position to disclose the physical meaning of Eq. (|5T7|). Taking up Eq. ( p.7[ ) 
with upper signs and setting p° on the energy-shell of the normal mode, p° = ^(X; p), we 
obtain, using the formulae in Appendix C, 
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(Z n + (X;P) 



-i 



+ 



<9Reo+(X;P) P) 



dP„ 



n(X-P) 

loop 



i a 



K 



(X;P). 



Here v" (= du™(X; p)/9p) is the group velocity of the mode [cf. Eq. (|C.2| )1. As will be 



loop 



shown in Appendix D, i \ &k) on the RHS is related to the net production rate, r™ et p , of 
the mode p° = ou™(X] p). Using Eqs. ( p.l| ) and ( |CJ.3| ) , we obtain 



9 n , v \ ® 

^ + v» + (A-;p).- 



+ 



N(X;P) 
dul(X;j>) 8N(X;P) 



dX» 



dP„ 



r net»(^; P) 



This can further be rewritten in the form, 



<V?pf;p) dn 



r netp(^; P) 



Similarly, Eq. 



dX dP 

(with upper signs) with p = — uj\ (X; — p) yields 
_ duj a + (X;p) dn _- a 

dx dp L »*p^>v)- 



Equation (|5.7|) with lower signs yields, on the energy-shell p = —uj t 1(X; — p), 



&j"(X;p) dn 
&X dP 

and on the energy-shell po = tot (X; p), 



net p 



(X;p) 



^ + v^(.Y;p).A)„(A-;^(A-;p),p) 



du a _(X]p) dn 



r netp(^; P)- 



(5.8) 



(5.9) 



(5.10) 



(5.11) 



<9X dP 

Equations ( p.8| ) - (|5.11 ) are the generalized relativistic Boltzmann equation for (anti)fermi- 
onic quasiparticles. 
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VI. PERTURBATION THEORY 



As has been discussed in the preceding section, the propagator in the physical- N scheme 
is free from the pinch singular term (in the narrow-width approximation) and then the 
perturbative calculation of some quantity yields "healthy" perturbative corrections. For 
constructing a concrete perturbative scheme, one more step is necessary. 

To extend the condition ( |5.7| ) to off the energy-shell, we divide / into two pieces [cf. 
Eq. Q] 

f(X; P) = 9(-p°) + e(p°)N(X; P) 
= f(X;P) + f (X;P) 

= [0(-p°) + e(p°)N(X; P)) + e(p°)N (X; P). 

(6.1) 

fo (and then also /) is defined as follows: Let Hi(X; p) (i = 1, 2, ...) be a region in a p°-plane 
that includes ith energy-shell. We choose lZi(X; p), such that, for i ^ j, IZi fl TZj = 0. On 
each energy-shell, N (X; P) = N(X; P), and, in whole p°-region but Hi(X] p) (i = 1, 2, ...), 
A^ (X; P) vanishes. dN (X; P)/dX and dN (X; P)/dP exist and N (X; P) obeys 

P ± ■ d x N T i (4) l0 ° P ± {Rea± JV } p B = 0. (6.2) 

Then, GP is given by Eq. ( |4.13|) with / for /. It is obvious from the above construction 
that this G$ does not possess pinch singularities in narrow-width approximation, and thus 
"healthy" perturbation theory is established. 

It is worth mentioning that there is arbitrariness in the choice of the regions , IZ^X; p) 
(i = 1, 2, ...). Furthermore, the choice of the functional forms of / and of fo is also arbitrary, 
provided that 

f(X° = Xl X; P) + f (X° = Xl, X; P) 
= /(X° = At,X;P), 
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where f(X° = X® n , X; P) is the initial data with Xf n the initial time. As has been discussed 
at the end of Sec. Ill, these arbitrariness are not the matter. 

To summarize, to the gradient approximation the (resummed) propagator Q of the theory 

is 

g = G + 5i + S 2 + S 3 . (6.3) 



Here Si - S 3 are as in Eqs. ( |2.7|) , (|A.1| ) - ( |A.2| ), and G is as in Eq. ( |4.S| ) provided that G 



(0) 
K 



is given by Eq. ( |4.13|) with / for /. / consists of two pieces as in Eq. ( |6.1| ). f (= e(p°)N ) 
subjects to Eq. ( |6.2|) , which is to be solved under a given initial data. It is to be noted that 
Si - S3 are the nonleading parts, so that, if one wants, for / in Si - S3, one can substitute 
fs, the solution to the "free Boltzmann equation," Eq. ( |2.13| ). 

As in m , determination of / or N proceeds order by order in perturbation theory, which 
we do not reproduce here. 

The vertex factor and the initial correlations are the same as in standard CTP formalism, 
except that an additional two-point vertex iL c (x,y)A_, Eq. ( ^6|) , exists. It is to be noted 
that the two-point vertex iL c (x,y)A_ has been built into G (cf. Eqs. fl4.12| ) and fl4.13j ) with 
Eqs. (|B.9|) and ( B.7|) ) and is absent in the perturbative framework using Q in Eq. ( |6.3|) . 



VII. SUMMARY AND DISCUSSION 

In this paper we have dealt with out-of-equilibrium perturbation theory for massless 
Dirac fermions. The fermion propagator is constructed from first principles. Essentially 
only approximation we have employed is the so-called gradient approximation, so that the 
perturbative framework applies to the quasiuniform systems near equilibrium or the nonequi- 
librium quasistationary systems. The framework allows us to compute any reaction rates 
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There comes out naturally the generalized Boltzmann equation (GBE) that describes the 
spacetime evolution of the number densities of quasiparticles, through which the evolution 
of the system is described. 
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Comparison of the present formalism with some earlier works has been made in |J . We 
like to add here two related works |I^,|T^], in which fermions are dealt with. In these papers, 



the GBE is derived in a traditional manner by starting with the SD equation (cf. Eq. ( |4.6| )). 
As is mentioned in |3|], the SD equation is nothing more than an equation that serves as 
resumming the self-energy part to makeup the resummed propagator, Eq. ( |4.8| ). Then, in 
order to derive the GBE, an additional input or condition is necessary. Our condition is 
Eq. ( |5.7|) . The additional input in fl3| , |l4|1 is, in our notation, an introduction of /' through 



Gi 2( 2i)(X; P) ~ 2mf'(X; P)e{p )fa{P) 
= 2m [e( Po )N'(X;P) 

-e{- VQ )(l-N'{X-P))]ta{P) 

(7.1) 

with cr(P) — 8(P 2 ). In fl3f , this is done on the basis of the quasiparticle picture, and, 
in [14j], on the order-of-magnitude estimation of E in diagrammatic analysis. Substituting 
Eq. ( |7. 1| ) into the Kadanoff-Baym equation, which is a part of the SD equation, the GBE 
for N'(X; P) and N'(X; P) results. In contrast to the formalism in the present paper, no 
counter Lagrangian is explicitly introduced there and then the consistency check for the 
formalism seems to be necessary. 

In [0,|J , comparison has been made of the present formalism with nonequilibrium thermo 
field dynamics (NETFD) p|. (See also pj). We like to recapitulate here that, in NETFD, 



a counter action is introduced on the basis of the renormalizability argument. Imposition of 
the renormalization condition leads to the GBE. 

The "derivation" in this paper of the GBE for a quark-gluon plasma (nonequilibrium 
QCD) is quite different from the traditional derivation (cf. [[L3|,[14]] ) . What we have shown 
here is that the requirement of the absence of large contributions from the perturbative 
framework leads to the GBE. This means that the quasiparticles thus defined are the well- 
defined modes in the medium. Conversely, if we start with defining the quasiparticles such 
that their number density functions subject to the GBE, then, on the basis of them, well- 
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defined perturbation theory may be constructed. In |12], the GBE is derived by imposing 
the renormalization condition" for the propagator. Since the pinch singularity, which arises 
in the narrow-width limit in our formulation, is a singularity in momentum space, it is 
not immediately obvious how to translate this condition into (space)time representation, as 
adopted in jl^]. Nevertheless, closer inspection of the structures of our formalism and of the 
NETFD tells us that our condition is in accord with on-shell renormalization condition in 
NETFD. 
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APPENDIX A: EXPRESSION FOR THE NONLEADING PART OF S(X; P) 
The nonleading parts, S 2 and S3, of S read 



S 2 (X; P) = - [A R (P) - A A (P)} e{p°) lb tN { -\P)A^ 

'a^W^p 



(A.l) 



S 3 (X;P) = -[A R (P)-A A (P)}e(p° 
+ LV -e 2ij a 0i p 



p°+p 3 j v 1 



p° + p 3 



(A.2) 



where e 123 = 1, = ~ 7 / V)/2, and 

N(-\X;P) = i [9(p°) {N ++ (X; \p°\,p) - AL_(X; \p°\,p)} 

-6(-p°) {iV ++ (X; \p°\, -p) - iV„_(X; \p% -p)} 
~ N { -\P), 

at(p) = e(p°)N„ + (\p \,p) + e(-p°)N- + (\ Po \, -p). 



As mentioned at the beginning of Sec. IIB, we have assumed that S 2 and S3 are, at most, 
of the same order of magnitude as the gradient term S\, and then X-dependence of N^~' 
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and Af have been ignored. It can readily be seen that, to the gradient approximation, 
if 3 x S2{x,y) = if 3 x S 3 (x,y) = holds. One can see from Eq. ( |A.2j ) that S 3 does not 
contribute to the unpolarized quantities. 

For a charge-conjugation- invariant system, N^~'(P) = N^(—P) and N(P) = J\f(—P) 
hold. 



APPENDIX B: SOLVING THE SD EQUATION, EQ. ( |4.6| 



Multiplying B L l - = t 3 B R f 3 - (-Br 1 = -f 3 B L f 3 ) [cf. Eq. (fO|] from the left (right) of each 
term in Eq. (|4.6| ), we obtain 



G = S + S-Z-G = S + G-Z-S, 



G = B L 1 -G- B R , 



S = B L 1 -S- B R l 



( 



Tj = B& ■ S ■ Bi 



ifr^R S K 
-i$A A 
\ 



V 



-Z A 

Here T, R ^ A ) = £n + £12(21), an d Sk and E# are the inverse Wigner transforms of 



/ 



S K (X;P) 



} xf(X;P) 



7° + 7 3 

pO _|_ p3 



P-d x f{X;P) 



P 

P2- 



and of 



S^(X;P) = [1-/(X;P)]S 12 (X;P) 
+/(X;P)E 21 (X;P) 



(B.l) 
(B.2) 

(B.3) 
(B.4) 



(B.5) 



+ 2 + /}rb. > 



(B.6) 
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respectively. Here {..., ...} PB is as in Eq. (|L4j). In obtaining Eq. ( |B.5|) [Eq. ( |B.6|) 1, use has 
been made of Eq. ( [4.5| ) [Eq. (|4.7|) 1. Although the last term in Eq. ( |B.6| ) may be dropped to 
the approximation under consideration, we have kept it. 

From the above definitions of £#(.4), ^k, and of (Eq. ([31])), we have 



v( c ) — v( c ) — n 

^R — — u > 

^K ) (X;P) = l ^±^P.d x f(X;P). 

pU _|_ 



(B.7) 



As seen from Eqs. ( |B.3| ) and 
may easily be solved to yield 

/ 

G = 



S and £ are triangular matrices, so that Eq. (|B.1|) 



^R 



V 



o 



-G A 



G 



R(A) 



G' K — Gr 



(i0A R ) '-Sk-^Aa) X -Y. k 



G, 



(B.8) 



(B.9) 



Substituting this back into G = Bl-G-Br, Eq. (P-2| ), we obtain, after Wigner transformation, 
Eq. ((4l| with 



G K (X; P) = G' K {X- P) + - {Gr + G A , f} PB . 



From Eq. 



, we obtain, after some manipulation, Eq. ( |4.9p with Eq. ( |4 . 1 U| ) and with 

G ( rU(X; P) = i [g+ {A) (X; P) - g R{A) (X; P) 

x E [757 ■ (p x Vx(r R(A) (X; P)' 



T=± 



X 



9r { a)(X;P) 



(B.10) 



where g R{A) is as in Eq. (f4.11|) and a T R{A) = a\ x + a[ 2(21) . 
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APPENDIX C: ON THE ENERGY SHELL 



Here we display some formulae, which hold on the energy-shells of quasiparticles [cf. 
Eqs. (|5.1|) and ( |5.2j ) with Eqs. (|4.10|) and ( 4.11f) 1 . In most formulae in this Appendix, we 
drop the argument A. 
Normal modes: 

We define the wave-function renormalization factors through taking derivative of 
Eq. ( |5.1| ) with respect to p°: 



(z£(t4(p) )P )) 1 



i=f 



dRea^(p°,±p) 



dp 



The group velocities of the modes are obtained from the definition (|5.1| ), 



v±( P ) 



dp 

2£04(p),p) 

dRea±(p°,±p) 



P + 



dp 



p° = ±Wj(p) 



By differentiating Eq. ( |5.1| ) with respect to A, we obtain 

<M(A;p) _ «9Rea+(A;^(A;p),p) 



<9A 



<9A 



(C.l) 



(C.2) 



(C.3) 



Abnormal modes: 



(z»(u4(p), P ; 



-i «9Rea|(p°,Tp) 



IT 



9 °= T u4( P ) 



v±(pj 



-^K( P ), P ) 

dRea±(p°,Tp) 



P + 



<9p 



p° = =F^(p) 



(C.4) 
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APPENDIX D: NET PRODUCTION RATES 



p° = u± a (p) fermionic modes: 



From Eqs. ( 4.10 ), (5.1), and ( |5.2j) , we see that the projection operators onto p° = lu'^ 
(p° = uot) mode is lf > + (#-)■ Then, the production and decay rates are written, in respective 
order, as 171,01 



K /a (p) 



* rvnla, n/a, \ \ 

-Z±' (w ± ' (p),p) 



xTr 



Si 2 (^ /a (p), P ))'^/ ± 



loop 



T n J a (X;p) 



±iZl /a (of 2 (w± /a (p),p 

^ r^n/a, n/a,\ \ 
~^ Z ± {V± (p),P) 



xTr 



loop 



(s 21 (^"(p),p))"" ,p / 



TiZl /a (a^(w± /a (p),p; 



loop 



where Z's are the wave-function renormalization factor, Eqs. ( |C.1| ) and ( p.4j ). Thus, the net 
production rate is 



where P * s as i* 1 Eq- (|4.15|) . 

p) antifermionic modes: 



net p' 



p) = [l-n(^ /a (p),p)]r^(X;p 
-n(^ /a (p),p)C/ a (p) 
= ±zi^>^(p),p) 



x(4(^ /a (X;p), P ; 



loop 



n /a 
P = — Cg> t ' 



(D.l) 



r » (-P) 



* rvnla, n/a, \ \ 

--Z T < [uj^ (-p),-p) 



xTr 



(E 21 (-^(-p),p)) too ^ 



loop 



T d (-P) 



^ rynla, n/a, \ \ 

-Z T < (u T < (-p),-p) 
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xTr 



(E 12 (-^(-p),p))""^ 



\ loop 



r n t P (-p) = [i - nK /a (- P ), - P )]r;-(-p) 



n/a i 



-n/a . 



-n( W f(-p),-p)rr(-p) 

±i^/>^ a (-p),-p) 

x(^(-^(-p),p)) lo ° P 



n/a , 
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